Bloch vectors for qudits 

Reinhold A. BertlmanEB and Philipp Krammeilll 

Faculty of Physics, University of Vienna, 
Boltzmanngasse 5, A-1090 Vienna, Austria 

We present three different matrix bases that can be used to decompose density 
matrices of d-dimensional quantum systems, so-called qudits: the generalized Gell- 
" Mann matrix basis, the polarization operator basis, and the Weyl operator basis. 

' Such a decomposition can be identified with a vector — the Bloch vector, i.e. a 

■ generalization of the well known qubit case — and is a convenient expression for 

^ . comparison with measurable quantities and for explicit calculations avoiding the 

^ I handling of large matrices. We present a new method to decompose density matrices 

via so-called standard matrices, consider the important case of an isotropic two- 
qudit state and decompose it according to each basis. In case of qutrits we show a 
representation of an entanglement witness in terms of expectation values of spin 1 
measurements, which is appropriate for an experimental realization. 



O^. I. INTRODUCTION 



The state of a (i-dimensional quantum system — a qudit — is usually described by a 
d X d density matrix. For high dimensions, where the matrices become large (for composite 
! systems of n particles the matrices are of even much larger dimension c?" x rf"), a simple 
way to express density matrices is of great interest. 

Since the space of matrices is a vector space, there exist bases of matrices which can be 
■ used to decompose any matrix. For qubits such a basis contains the three Pauli matrices, 
QO ! accordingly, a density matrix can be expressed by a 3-dimensional vector, the Bloch vector, 
and any such vector has to lie within the so-called Bloch ball [H, 0]- Unique for qubits is the 
fact that any point on the sphere, Bloch sphere, and inside the ball corresponds to a physical 
^ . state, i.e. a density matrix. The pure states lie on the sphere and the mixed ones inside. 

In higher dimensions there exist different matrix bases that can be used to express qudits 
as [d"^ — l)-dimensional vectors as well. Different to the qubit case, however, is that the 
map induced is not bijective: not every point on the "Bloch sphere" in dimensions d"^ — 1 
corresponds to a physical state. Nevertheless the vectors are often also called "Bloch vectors" 
(see in this context, e.g., Refs. 0, 0, H, M, 01)- 

In this paper we want to present and compare three different matrix bases for a Bloch 
vector decomposition of qudits. In Sec. [TTl we propose the properties of any matrix basis for 
using it as a "practical" decomposition of density matrices and recall the general notation 
of Bloch vectors. In Sees. IIIII - |V] we offer three different matrix bases: the generalized 
Gell-Mann matrix basis, the polarization operator basis, and the Weyl operator basis. For 
all these bases we give examples in the dimensions of our interest and present the different 
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Bloch vector decompositions of an arbitrary density matrix in the standard matrix notation. 
Next in Sec. IVll by constructing tensor products of states we study the isotropic two-qudit 
state and present the results for the three matrix decompositions, i.e. for the three different 
Bloch vectors. In Sec. IVIII we focus on the isotropic two-qudit state and calculate the 
Hilbert-Schmidt measure of entanglement (see, e.g., Refs. 0,1, [13, [ill). Its connection to 
the optimal entanglement witness is shown, which is determined in terms of the three matrix 
bases. An example for the experimental realization of an entanglement witness is given in 
Sec. IVII B[ The mathematical and physical advantages/disadvantages by using the three 
different matrix bases are discussed in Sec. IVIIIt where also the final conclusions are drawn. 



II. PRELIMINARIES 

A qudit state is represented by a density operator in the Hilbert-Schmidt space acting 
on the d-dimensional Hilbert space Ti'^ that can be written as a matrix — the density 
matrix — in the standard basis {\k)} , with k = 1,2, . . . d or k = 0,1,2, . . . d — 1. 

Properties of a "practical" matrix basis. For practical reasons the general properties of a 
matrix basis which is used for the Bloch vector decomposition of qudits are the following: 

i) The basis includes the identity matrix 1 and d — 1 matrices {Ai} of dimension d x d 
which are traceless, i.e. TrAj = . 

ii) The matrices of any basis {A^} are orthogonal, i.e. 

Ti AlAj = N6ij with NeR. (1) 

Bloch vector expansion of a density matrix. Since any matrix in the Hilbert-Schmidt space 
of dimension d can be decomposed with a matrix basis {A^}, we can of course decompose a 
qudit density matrix as well and get the Bloch vector expansion of the density matrix, 

P=^l + &-f, (2) 

where 6 ■ F is a linear combination of all matrices {Ai} and the vector b G H'^^^^ with 
bi = (Fj) = TrpFj is called Bloch vector. The term ^1 is fixed because of condition Trp = 1. 

Remark. Note that a given density matrix p can always be decomposed into a Bloch 
vector, but not any vector a that is of the form ([2]) is automatically a density matrix, even 
if it satisfies the conditions Trcj = 1 and Tro"^ < 1 since generally it does not imply a > 0. 

Each different matrix basis induces a different Bloch vector lying within a Bloch hyper- 
sphere where, however, not every point of the hypersphere corresponds to a physical state 
(with p > 0); these points are excluded (holes). The geometric character of the Bloch space 
in higher dimensions turns out to be quite complicated and is still of great interest (see 

Refs. SHBBH). 

All different Bloch hyperballs are isomorphic since they correspond to the same den- 
sity matrix p. The interesting question is which Bloch hyperball — which matrix basis — is 
optimal for a specific purpose, like the calculation of the entanglement degree or the deter- 
mination of the geometry of the Hilbert space or the comparison with measurable quantities. 
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III. THE GENERALIZED GELL-MANN MATRIX BASIS 

A. Definition and example 

The generalized Gell-Mann matrices (GGM) are higher-dimensional extensions of the 
Pauli matrices (for qubits) and the Gell-Mann matrices (for qutrits), they are the standard 
SU(N) generators (in our case N = d). They are defined as three different types of matrices 
and for simplicity we use here the operator notation; then the density matrices follow by 
simply writing the operators in the standard basis (see, e.g. Refs. (sl. [l2|): 

d{d-l) 



i) ^ „ ' symmetric GGM 



Af = \j){k\ + \k){j\, l<j<k<d, (3) 

ii) ^^iftli antisymmetric GGM 

Af = -^\J){k\ +t\k){j\, l<j<k<d, (4) 

iii) {d - 1) diagonal GGM 



A' = Wtt^^ >^b)01 - + + h l<l<d-l. (5) 




In total we have d"^ — 1 GGM; it follows from the definitions that all GGM are Hermitian 
and traceless. They are orthogonal and form a basis, the generalized Gell-Mann matrix 
basis (GGB). A proof for the orthogonality of GGB we present in the Appendix I A 1[ 

Examples. Let us recall the case of dimension 3, the 8 Gell-Mann matrices (for a repre- 



sentation see, e.g., Refs. [111. Il 



i) 3 symmetric Gell-Mann matrices 

.12 _ I 1 n n 1 AP = I 1 , \f 






(6) 



ii) 3 antisymmetric Gell-Mann matrices 
/ -i 0\ 
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(7) 

ii) 2 diagonal Gell-Mann matrices 

= ( J -1 I , = I 1 I . (8) 
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To see how they generahze for higher dimensions we show the case we need for qudits of 
dimension d = 4: 
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ii) 6 antisymmetric GGM 
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iii) 3 diagonal GGM 
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Using the GGB we obtain, in general, the following Bloch vector expansion of a density 
matrix: 

p = ^1 + 6-A, (12) 

with the Bloch vector b = {{bi''}, {b{''}, {&'}) , where the components are ordered and for the 
indices we have the restrictions 1 < j < k < d and 1 < / < d — 1 . The components are given 
by bi^ = TrAi'^p , fe^'^ = TrA^'^p and 6' = TrA'p . All Bloch vectors lie within a hypersphere 
of radius \b\ < sj^d — l)/2d. For example, for qutrits the Bloch vector components are 
b = {bl'^,bl^,bf,bl^,bf,bf,b\b'^) corresponding to the Gell-Mann matrices (E]), (CD, 1^ 
and |6| < ^/ijs . 

As already mentioned the allowed range of b is restricted. It has an interesting geo- 
metric structure which has been calculated analytically for the case of qutrits by studying 

2- dimensional planes in the 8-dimensional Bloch space or numerically by considering 

3- dimensional cross-sections In any case, pure states lie on the surface and the mixed 
ones inside. 
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B. Standard matrix basis expansion by GGB 

The standard matrices are simply the dx d matrices that have only one entry 1 and the 
other entries and form an orthonormal basis of the Hilbert-Schmidt space. We write these 
matrices shortly as operators 

\3){k\, with j,fc = l,...,d. (13) 

Any matrix can easily be decomposed into a "vector" via a certain linear combination of 
the matrices f|T3l) . Knowing the expansion of matrices f|T3|l into GGB we can therefore find 
the decomposition of any matrix in terms of the GGB. 

We find the following expansion of standard matrices ( fT3l) into GGB : 



(14) 



' \ [h^J" + ikf) for j <k 

= j |(AJ^- - 2A^^) forj>fc 

A^-i + V ^ A-'+" + i 1 hi ] = k. 

2j ^2(j+n)(i+n+l) d J 

Proof. The first two cases can be easily verified. 

To show the last case we first set up a recurrence relation for |/)(/|, which we obtain by 
eliminating the term X]j=i ^"^"^ expressions ([5]) for A' and A'~^ 



and we consider the case I + 1 = d 



\d-l){d-l\ = -xI^tA-^^'-" + xhr^,^'-' + \d){d\. (16) 



2{d-l) \2{d-l) 



From A*^ ^ given by Eq. 



' d-l 



we get the Bloch vector decomposition of \d){d\ 



U .l{d-l)d^,_. 



where we have applied Yl'jJi b)OI = 1 — \d){d\ . 



Inserting now decomposition (1181) into relation ( 1161) we gain the Bloch vector expansion 
for \d — l){d — 1| and recurrence relation (|T5l) provides \d — 2){d — 2| and so forth. Thus 
finally we find 



\d-n){d-n\ = -x T , ^ Xd-n-i _^Sr ^ \<i-n+k^}^^ 

]l ^{d-n) ^ ^2{d-n + k + l){d-n + k) d' 

(19) 

the relation we had to prove, where d — n = j . □ 
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IV. THE POLARIZATION OPERATOR BASIS 

A. Definition and examples 

The polarization operators in the Hilbert-Schmidt space of dimension d are defined as 
the following d x d matrices [3, 14| : 



Tlm = J2 ^sZ LM I ^) I ■ (20) 



k,l=l 

The used indices have the properties 

„ _ cbi 
2 ' 

L = 0,1, ... ,2s, 
M = -L, -L + 1,...,L-1,L, 
mi = s, m2 = s — 1, . . . , rrid = —s . (21) 

The coefficients C^™* are identified with the usual Clebsch-Gordan coefficients Cj^'L^^j^^ 



of the angular momentum theory and are displayed explicitly in tables, e.g., in Ref. [14 



For L = M = the polarization operator is proportional to the identity matrix 0, 14 

Too = 1 . (22) 

It is shown in Ref. |J] that all polarization operators (except Too) are traceless, in general 
not Hermitian, and that orthogonality relation ([1]) is satisfied 

'^^'^LiMi'^L2M2 = ^LiL2^MiM2 ■ (23) 

Therefore the d"^ polarization operators fl20l) form an orthonormal matrix basis — the 
polarization operator basis (FOB) — of the Hilbert-Schmidt space of dimension d. 

Examples. The simplest example is of dimension 2, the qubit. For a qubit the FOB is 
given by the following matrices (s = 1/2; L = 0, 1; M = —1, 0, 1) 

^ _ 1 / 1 0\ „ fO 1 



For the next higher dimension d = 3 {s = 1), the case of qutrits, we get 9 polarization 
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operators T^a/ with L = 0, 1, 2 and M = —L, L and we have 
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r2„2 = I I . (25) 



Then the decomposition of any density matrix into a Bloch vector by using the POB has, 
in general, the following form: 

P = -J^ + J2 Yl ^lmTlm = ^1 + 6-f , (26) 

L=l M=-L 

with the Bloch vector h = &10, &11, ^2-2, ^lm), where the components are 

ordered and given by Blm = '^'^T'ImP ■ general the components Iilm are complex since the 
polarization operators Tlm are not Hermitian. All Bloch vectors lie within a hypersphere of 
radius \b\ < yj{d- l)/d. 

In 2 dimensions the Bloch vector h = (&1-1, &10, &11) is limited by \b\ < and forms a 
spheroid the pure states occupy the surface and the mixed ones lie in the volume. This 
decomposition is fully equivalent to the standard description of Bloch vectors with Pauli 
matrices. 

In higher dimensions, however, the structure of the allowed range of b (due to the posi- 
tivity requirement p > 0) is quite complicated, as can be seen already for = 3 (for details 
see Ref. j^). Nevertheless, pure states are on the surface, mixed ones lie within the volume 

— # — * 

and the maximal mixed one corresponds to |6| = 0, thus \b\ is a kind of measure for the 
mixedness of a quantum state. 



B. Standard matrix basis expansion by POB 

The standard matrices f|T3l) can be expanded by the POB as [li 



1^) 0' I = E E V 1^ LM Tlm . (27) 

L M ^ 



Note that YIim actually fixed by the condition rrij + M 



nii. 
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Proof. Inserting definition (l20l) on the right-hand side (RHS) of equation (1271) we find 



RHS 



j ^ '^L + l ^ ^^^^ 



kl 



2s + 1 



^ srrij , LM smi, LAI 



\k){l\ 



k,l 



where we used the sum rule for Clebsch-Gordan coefficients [14 

2c + 1 



21) -{- 1 '^'^ C7 



c,7 



(28) 



(29) 



V. WEYL OPERATOR BASIS 

A. Definition and example 

Finally we want to discuss a basis of the Hilbert-Schmidt space of dimension d that 
consists of the following (P operators: 

d-l 

Unm = ^e^'^"|A;)((A; + m)modrf| n,m = 0,1, ... ,d - 1 , (30) 

where we use the standard basis of the Hilbert space. 

The operators in notation fl30|) have been introduced in the context of quantum tele- 



portation of audit states [15[ and are often called Weyl operators in the literature (see e.g. 
Refs. 16|, llJ, ll8|). The d'^ operators fl30l) are unitary and form an orthonormal basis of the 
Hilbert-Schmidt space 

TlUljJij = d5nl5mj (31) 

(a proof is presented in Appendix lA 31) - the Weyl operator basis (WOB). They can be used 



to create a basis of rf^ maximally entangled qudit states 16|, ll9|, |20 |. 



Clearly the operator Uqq represents the identity Uoo = 1 

Example. Let us show the example of dimension 3, the qutrit case. There the Weyl 
operators ( 130|) have the following matrix form 



f/oi = 1 , Uo2 = \ l 0] , (32) 





'10 \ / 1 

Uw = \ e2-/3 , f/n = e^-/^ | ^ 
^0 e-^^'l^l 

'10 \ / 1 

f/20 = I e-2-/3 , f/21 = e-2-/3 I ^ 

g Q g2W3 / 
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Using the WOB we can decompose quite generally any density matrix into a Blocli vector 

1 '^"^ 1 

P = ^ 1 + ^nmUnm = ^1 + b-U, (33) 

n,m=0 

with n,m = 0,1, ...,d — 1 (600 = 0). The components of the Bloch vector b = ({&„m}) 
are ordered and given by bnm = Tr Unm p ■ In general the components bnm are complex 
since the Weyl operators are not Hermitian and the complex conjugates fulfil the relation 
Kim — ^~^^^b^n-m , which follows easily from definition (!30!) together with the hermiticity 
of p . 

All Bloch vectors lie within a hypersphere of radius |6| < ^/ d — 1/d . For example, for 
qutrits the Bloch vector is expressed by 6 = (6017 &02) ^lOi bu, 612, &20) ^2i) ^22 ) and |6| < V2/3. 
In 3 and higher dimensions the allowed range of the Bloch vector is quite restricted within 
the hypersphere and the detailed structure is not known yet. 

Note that in 2 dimensions the WOB as well as the GGB coincides with the Pauli matrix 
basis and the FOB represents a rotated Fauli basis (where <7± = ^ (cxi ± 1(12)), in particular 



{f^oo,f^oi,f^io,f^ii} = {l,o-i,o-3,«cr2} , (34) 
{l,Af,Af,Ai} = {l,a,,a2,as} , (35) 

{Too,Tii,Tio,Ti_i} = j-^ 1, -0-+, -^0-3, . (36) 

B. Standard matrix basis expansion by WOB 

The standard matrices ( fT3l) can be expressed by the WOB in the following way 

d-l 

= ;^ (37) 



Proof. We insert the definition of the Weyl operators ( l30i) on the right-hand side (RHS) 
of Eq. (ETD, use Eq. (IAM and get 

^ d-l 

RHS = - J] e^'('-^) |r)((r + A;-j)niodrf| 

l,r=0 

d-l d~l 

= \j){k\ + Yl $^e^^('-^')|r)((r + A;-j)niodd| 

r^j, r=0 1=0 

= □ (38) 

VI. ISOTROPIC TWO QUDIT STATE 

Now we consider bipartite systems in a. d x d dimensional Hilbert space Ti.^ The 
observables acting in the subsystems Ha and Hb are usually called Alice and Bob in quantum 
communication. 
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Quite generally, a density matrix of a two-qudit state acting on ®7iJ can be decom- 
posed in the following way (neglecting the reference to A and B) 



p = - 1 (g) 1 + ni (g) 1 + mi 1 (g) Ti + (g) r 

a 



Hi, rrii, Cij G C 



(39) 



where {Fj} represents some basis in the subspace T-C^ . The term Cij Fj (g F^ always can be 
diagonalized by two independent orthogonal transformations on Fj and Tj [21]. Altogether 
there are (d^)^ — 1 terms. 

However, for isotropic two-qudit states — the case we consider in our paper — the second 
and third term in expression (!39l) vanish and the fourth term reduces to cu F j ® F j , which 
implies the vanishing of (c/^ — l)^ + ((i^ — 1) = (i^((i^ — 1) terms. Consequently, for an 
isotropic two-qudit density matrix there remain d"^ — 1 independent terms, which provides 
the dimension of the corresponding Bloch vector. Thus the isotropic two-qudit Bloch vector 
is of the same dimension — lives in the same subspace — as the one-qudit vector, which is a 
comfortable simplification. 



Explicitly, the isotropic two-qudit state p^a'' is defined as follows [22|, |23|, [24 



a 



rf2 



1 , a eR. 



d^ - 1 



< a < 1 



where the range of a is determined by the positivity of the state. The state 
state, is maximally entangled and given by 



(40) 
6^), a Bell 

(41) 



where denotes the standard basis of the d-dimensional Hilbert space. 



A. Expansion into GGB 

Let us first calculate the Bloch vector notation for the Bell state |(/)^|_) {(f)'^\ in the GGB. 
It is convenient to split the state into two parts 



j,k=i 

A + B, 



d 



where A and B are defined by 



j<k 



j<k 



(42) 

(43) 
(44) 



and to calculate the two terms separately. 
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For term A we use the standard matrix expansion (fT^ for the case j ^ k and get 



.j<k 



j<k 



1 J2 (Af ® Af - Af ® Af ) . 



(45) 



For term B we need the case j = A; in expansion ( |T^ and obtain after some calculations 
(the details are presented in Appendix IA2p 



d-l 



B 



l5^A™®A- + ll®l 

m=l 



(46) 



Thus all together we find the following GGB Bloch vector notations, for the Bell state fl42|) 

and for the isotropic two-qudit state (HOjl 

. i_l«l + iiA. (48) 

where we defined 

d-l 

A := ^ Af ® Af - Af ® Af + ^ A" ® A'" . (49) 



i<j 



i<j 



m=l 



B. Expansion into POB 

Now we calculate the Bell state \4>'\.) in the POB. Using expansion ( |27l) and the sum 
rule for the Clebsch-Gordan coefficients 14 



we obtain 



9r 4- 1 



a,7 



26 + 1 



(50) 



1 



1 ^ v/(2L + l)(2L- + l) / ^ \ 

^ 2s + 1 I ^ ^smj,LAf'-"smj,L'M I 

L,L' ' \ ij / 

l^ v/(2L + l)(2L- + l) 



2L + 1 

L,L' 



(51) 



12 



where we extracted the unity (recall Eq. (l22l) ) and defined 

T := Yl Tlm®Tlm. (52) 

L,M^O,0 

Result (ISTl) provides the POB Bloch vector notation of the isotropic two-qudit state ( l40l) 

pL'^) = ^1®1 + ^T. (53) 

C. Expansion into WOB 

Finally we present the Bell state in the WOB (the details for our approach using the 
standard matrix expression (137|) can be found in the Appendix lA 41 see also Ref. [3]) 



|0t>(0+| = J2^®^ + J2U. (54) 



with 



U := Y,Uira®U^lra, (/, m) ^ (0, 0) , (55) 

/,m=0 

where negative values of the index / have to be considered as mod d , and from formula (15^ 
we find the WOB Bloch vector notation of the isotropic two-qudit state 

p(f) = + (56) 



VII. APPLICATIONS OF THE MATRIX BASES 
A. Entangled isotropic two-qudit states 



In Ref. [ll| the connection between the Hilbert-Schmidt (HS) measure of entanglement 
ITol | and the optimal entanglement witness is investigated. Explicit calculations for 
both quantities are presented in case of isotropic qutrit states. For higher dimensions, the 
isotropic two-qudit states, the above quantities are determined as well but in terms of a 
rather general matrix basis decomposition. With the results of the present paper we can 
calculate all quantities explicitly. Let us recall the basic notations we need. 

The HS measure is defined as the minimal HS distance of an entangled state pent to the 
set of separable states S 

D{pent) ■= min||p-pent|| = ||P0-Pent|| , (57) 

where po denotes the nearest separable state, the minimum of the HS distance. 

An entanglement witness A E A {A = Aa®Ab , the HS space of operators acting on the 
Hilbert space of states] is a Hermitian operator that "detects" the entanglement of a state 
Pent via inequalities [lO, 25, [2^, 27 . 
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Definition 1. An entanglement witness A is a Hermitian operator with the following proper- 
ties: The expectation value of A is negative for an entangled state, whereas it is non-negative 
for any separable state. 

{pent, A) = Trpent^ < 0, 

{p,A) = TTpA > WpeS. (58) 

The fact, however that there exists an operator satisfying inequahties fl58l) for any 
entangled state, i.e. that the definition is meaningful, has to be proved; it follows from 
the Hahn-Banach Theorem of functional analysis (for a simple geometric approach, see 



Ref. [11|) 



An entanglement witness is "optimal" , denoted by Aopt , if apart from Eq. (!58ll there 
exists a separable state po E S such that 

(po,Apt) = 0. (59) 

The operator Aopt defines a tangent plane to the set of separable states S and all states Pp 
with {pp, Aopt) = lie within that plane; see Fig. [TJ 

Let us call the lower one of the inequalities (158|1 an entanglement witness inequality, short 
EWI. It detects entanglement whereas a Bell inequality determines non-locality. Rewriting 
Eq. ([58D as 

{p,A) - (Pent, A) > \/peS, (60) 

the maximal violation of the EWI is defined by 



5(pent) = max min(p,A) - (pcnt,A) , (61) 

A, ||A-al||<l y pes y 

where the maximum is taken over all possible entanglement witnesses A, suitably normalized. 



Then an interesting connection between the HS measure and the concept of entanglement 
witnesses is given by the Bertlmann-Narnhofer-Thirring Theorem, illustrated in Fig. [1] [13] . 
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Theorem 1. 

i) The maximal violation of the EWI is equal to the minimal distance of pent to the set S 

B{pent) = D{p,^t) . (62) 

a) The maximal violation of the EWI is attained for an optimal entanglement witness 

^ ^ ^ Po - Pent - (P0,P0 - Pent) 1 ^g^^ 
IIPO — Pent II 

Thus the calculation of the optimal entanglement witness Aopt to a given entangled state 
Pent reduces to the determination of the nearest separable state po • In special cases po is 
detectable but in general its detection is quite a difficult task. We are able to find the nearest 
separable state by working with Lemma [H a method we call guess method [ll| . 

Lemma 1. A state p is equal to the nearest separable state po if and only if the operator 

^ ^ P - Pent - {P,P- Pent) 1 /g^x 
Hp - Pent II 

is an entanglement witness. 

Lemma^ [1] probes if a guess p is indeed correct for the nearest separable state. Then 
operator C represents the optimal entanglement witness Aopt (l63l) . 

Now let us apply the matrix bases we discussed in the previous sections and calculate 
the quantities introduced above. As an entangled state we consider the isotropic two-qudit 
state p^a^'^^^, that is the state pi'^'* (1401) for < a < 1. 

Starting with the GGB we can express that state in our Bloch vector notation by formula 
IHjl . By using Lemma [1] we find that the nearest separable state is reached at a = 

Pr-Pil^-;^1«1 + ^^A. (65) 



It provides the HS measure 



and the optimal entanglement witness (|63 



^op.(pi';L) = y^ISI-^^A, (67) 

where we used the HS norm ||A|| = 2y/W^. 

Clearly, the maximal violation B of the EWI equals the HS measure D 

B(n^'^^ ) - -in^'^^ A 

\Ha, ent) ~ \ Pa, ent ) ^opt 

Ad) 



d V d+1 



« - :rvT = D{pZnt) ■ (68) 
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For expressing above quantities by the matrix bases FOB and WOB it suffices to calculate 
the proportionality factors between A, T and U . By comparison of the three forms for the 
isotropic qudit state (l48l) . fl53|) and (!56|) we find 

A = 2T and T = . (69) 

(Jj 

It provides the following expressions, for the FOB 

+ (70) 



(d) . J- ld-1 



and for the WOB 



A.MZJ--,^—,tr^l-^=T, (71) 



^„„.(pSL) = (73) 

Of course, the HS measure -D(pQ'^ent) remains the same expression fl66l) independent of the 
chosen matrix basis, which can easily be verified using ||T|| = \/(P^-l and \\U\\ = d\/ — \ . 



B. Entanglement witness representation for experiments 

Entanglement witnesses are Hermitian operators and therefore observables that should 
be measurable in a given experimental set-up and thus provide an experimental verification 
of entanglement. The quantity to be measured is the expectation value 

(A) = TrAp (74) 

of an entanglement witness A for some state p. If (A) < then the state p is entangled. 
But which measurements have to be performed? 

Obviously it is appropriate to express the entanglement witness in terms of generalized 
Gell-Mann matrices (l3])-([5]), since they are Hermitian. For d = 3 — qutrits — the Gell-Mann 
matrices ©-([H]) can be expressed in terms of eight "physical" operators, the observables 



Sn,, Sy, S^, S^, Sy, {S^, Sy}, [S y , S'J, {S ^ , S of & splu-f systcm, where S = {S^, Sy, S^) is 
the spin operator and {Si, Sj} = SiSj + SjSi (with i,j = x, y, z) denotes the corresponding 
anticommutator. The decomposition of the Gell-Mann matrices into spin-1 operators is as 
follows (for a similar expansion, see Ref. 0]): 

= 21 + ^{hS.-3Sl-3S'y), X' = 1= (^-21 + ^{hS. + Si + S'y)y (75) 
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All operators can be represented by the following matrices: 




{S,,Sx} = ^ 1 -1 . (76) 



Thus we can express any observable on a n-qutrit Hilbert space — a composite system of n 
particles with 3 degrees of freedom — in terms of above spin operators (1761) . 

As an example we want to study the entanglement witness for the isotropic two-qutrit 
state, i.e. state (HOl) for d = 3. In this case we obtain for the optimal entanglement witness 

Aso = (l® 1 - ^A), (77) 



3^2 V 4 

(i.e. Eq. ( 1671) for d = 3) where the operator A is defined in Eq. (H9l) . 

Expressing the Gell-Mann matrices in A (H9l) by the spin operator decomposition (1751) we 
find for the expectation value of the entanglement witness Ajgo 



(A„> = ^{1S1>- j^(A), (78) 



where 



+ ^ {{Si ® si) + (5J ® s'^)) + ^ {{Si ® 5J) + (5J ® si) 

+ i {{{Sz, Sx} ® {S,, Sx}) - {{Sy, 5J ® {Sy, 5J) - {{Sx, Sy} ^ {Sx, Sy})^ . (79) 



Decomposition (1791) has to be determined experimentally by measuring the several expecta- 
tion values with the set-ups on both Alice's and Bob's side. 

The advantage of the entanglement witness procedure is that for an experimental outcome 
(Aiso) < the considered quantum state is definitely entangled, whereas in case of Bell 
inequalities a violation detects nonlocal states. That means by the entanglement witness 
procedure we are able to detect more entangled states than with Bell inequalities. The 
amount of measurement steps necessary to determine an entanglement witness is about the 
same as in the Bell inequality procedure (see, e.g., Refs. 28, 29, [s^, 31])- 
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VIII. CONCLUSION 

In this paper we present three different matrix bases which are quite useful to decompose 
density matrices for higher dimensional qudits. These are the generalized Gell-Mann matrix 
basis, the polarization operator basis and the Weyl operator basis. Each decomposition we 
identify with a vector, the so-called Bloch vector. 

Considering just one-particle states we observe the following features: The generalized 
Gell-Mann matrix basis is easy to construct, the matrices correspond to the standard SU(N) 
generators (A^ = d), but in general (in d dimensions) it is rather unpractical to work with 
the diagonal matrices due to their more complicated definition. On the other hand, the 
Bloch vector itself has real components, which is advantageous, they can be expressed as 
expectation values of measurable quantities. For example, in 3 dimensions the Gell-Mann 
matrices are Hermitian and the Bloch vector components can be expressed by expectation 
values of spin 1 operators. The polarization operator basis is also easy to set up, all you 
need to know are the Clebsch-Gordan coefficients which you find tabulated in the literature. 
However, the Bloch vector contains complex components. For the Weyl operator basis the 
corresponding operators are again simple to construct, they are non-Hermitian but unitary. 
The Bloch vector itself has a very simple structure, however, with complex components. Let 
us note that in 2 dimensions all bases are equivalent since they correspond to Pauli matrices 
or linear combinations thereof. 

In case of two-qudits we have studied the isotropic states explicitly and find the following: 
In the generalized Gell-Mann matrix basis the Bloch vector fH5]) with expression is more 
complicated to construct, in particular the diagonal part B fHBl) (see Appendix IA2p . In the 
polarization operator basis the Bloch vector flS5]) with expression (1321) can be easily set up 
by the knowledge of the Clebsch-Gordon coefficient sum rule fl5U]) and in the Weyl operator 
basis the Bloch vector (156|1 with definition (!55l) is actually most easily to construct. 

The Hilbert-Schmidt measure of entanglement can be calculated explicitly for all isotropic 
two-qudit states and we want to emphasize its interesting connection to the maximal viola- 
tion of the entanglement witness inequality. Theorem [H 

For the experimental realization of an entanglement witness the generalized Gell-Mann 
matrix basis is the appropriate one since the generalized Gell-Mann matrices are Hermitian. 
For a different task, however, the determination of the geometry of entanglement the Weyl 
operator basis turns out to be optimal. In our example of the entangled isotropic two- 
qutrit state the entanglement witness can be expressed by experimental quantities, the 
expectation values of spin-1 measurements. In this way one can experimentally find out 
whether a state is entangled or not, i.e., we can obtain rather precise information on the 
quality of entanglement. 

Quite generally, the Bloch vector decomposition into one of the three matrix bases is 
of particular advantage in the construction of entanglement witnesses. It turns out that 
if the coefficients of the decomposition satisfy a certain condition the considered operator 



represents an entanglement witness, i.e. satisfies inequahties flSSl) (for details see Ref. |32| ) 
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APPENDIX A 

1. Proof of Orthogonality of GGB 

We want to proof condition dH) for the GGB which consists of the d? - 1 GGM (El), (il), 
([5]) and the dx d unity 1. Since all GGM are Hermitian (thus TrA^Aj = TiA^Aj = TiAjAi) 
it suffices to proof the following conditions: 

(Al) 
(A2) 
(A3) 
(A4) 
(A5) 
(A6) 

Proof of condition (lAip . Inserting definition ([3]) we have 

d 

TrAfAr = J2(^\{\j){k\ + \k){j\){\m){n\ + \n){m\)\l) 
1=1 

= J]((/|j)(fc|m)(n|0 + {l\j){k\n){m\l) + m{j\m){n\l) + m{j\n){m\l)) 
I 

= S^^S''"^ -\- ^i^^^^ _|_ ^k'^^j^ _|_ ^km^jn 

= 2 5^''"5^", (A7) 
where we used in the last step that = since we have j < k and m < n. 

Proof of condition flA2p . This case is equivalent to the one before apart from changed 
signs that do not matter 



TrAf A^'^ 


= 2(5^"'(5*^" 


Tr A-^'^A"^" 

a a 


= 2(5-''™5^" 


TrA'A'" 


= 2(5'™ 


TrAf A^" 


= 


TrAf A'" 


= 


TrAf A"^ 


= 0. 



Tr A-''^A™" = — (5-'"5'^™ -|- ^i'^^k^ _|_ ^^^^j"^ _ ^km^jn 

jm ^kn 



= 2(5^^(5'^". (AS) 
Proof of condition (1A3I) . Using definition ([5]) and denoting 



where / < m without loss of generality, we get 

d I m 

TrA'A™ = CiC.m^(^^^{p\k){k\n){n\p) + lm{p\l + l){l + l\m + l){m + l\p) 

p=l k=l n=l 

I m 

-m^{p\k){k\m+l){m+l\p) - l^{p\l + + l\n){n\p) 

k=l n=l 

(I m \ 

/ + /m(5'" - m^(5'^(™+i) - j . (AlO) 

k=l n=l / 
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Using the fact that = for m > and 



= <'° '^^ = '^ (All) 

^ 1 Z if / < m ^ ^ 



n=l 

we obtain 



TrA^A"^ = + 1)5'"^ = 2 5'"^. (A12) 

Proof of condition ( 1A4I) . Analogously to the proofs ( 1A7I) and ( lASIl we find 

TrAf A™" = i (-^J"^*--'™ + - + ^i"^'^™) = 0. (A13) 

Proof of condition ( lASP . Inserting definitions ([3]) and gives 

d 

TrAfA™ = Cm'^(^-m{p\k){j\m + l){m + l\p) - m{p\j){k\m + l){m + l\p) 



p=i 

m 



+ '^{p\j){k\n){n\p) + ^{p\k){j\n){n\p) 



n=l n=l 

m 



1=1 

= 0, (A14) 
since per definition we have j < k . 

Proof of condition (\A6\i . This proof is equivalent to the previous one since constant 
factors in front of the terms do not matter. 

2. Calculation of term B in GGB 

To obtain the Bloch vector notation of term B (l44l) we insert the standard matrix ex- 
pansion (HM for the case j = k. We split the tensor products in the following way 

B = M fii + ^2 + ^3 + ^4 + i 1 ® 1 j , (A15) 
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where the terms _Bi, . . . , S4 are introduced by (note that A° = 0) 

. T d-j-1 
j=l V n(=0=0 



5i = V — -A^~' ® + V -A^+" ® A^+" (A16) 

lf\ 2j A. 2(j+n)(i + n + l) ^ ^ 



^2 = y I -y\l—, — — -A^'-i®A^-+' 



d-j-i 



J -I 



^ y 4j(j+n)(j+n + l) 

+ V — ® A^'+' (A17) 

n^W=o V(J +n){j+n + + + / + 1) ' 

^3 = ly" ( -j:^A^-i® 1 + y , ^ A^'+" 1 ) (AI8) 



^4 = -,y(-\- — -l^A'-' + V , ^ 1 ^ A^'+' 1 . (A19) 



d 

Only the first term Bi (IAT6D Gfives a contribution 



Si = ( ^ ^ + ^ ) A™ ® A™ = - V A"^ ® , (A20) 

^ \2(m + 1) 2m(m + 1) J 2^ ' ^ ^ 

m=l ^ ^ ' ^ ' ' m=l 



^ \2(m + 1) 2m(m + 1) 
whereas the remaining terms vanish: 

d-l 



m<p, m,p=l 



B2 = y -J-, — — , + , ^ A-^A^ 



m>p, m,p=l 



+ y I-a/tT TT^T TT+ , ^ )a-0A^ 



4(p + l)m(m+l) ^4p(p + l)m(m + 1) 



> — — + y — — A™ (g) A*' 

^m<p 2Vm(m + l)p{p + 1) 2v/m(m + l)p{p + 1) J 

= 0, (A21) 
and in quite the same manner 

, d-l 

53 = - V , A™ ® 1 = , 

^1 V2^(m + 1) 



= -y , 1 ® AP = . (A22) 



Thus we find the following Bloch vector of B (1441) 



1 "^"^ 1 
B = _^A™®A™ + -1®1. (A23) 

m=l 
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3. Proof of Orthonormality of WOB 

For proofs relevant in the WOB we often need the equivalence 

n=0 k 

So we use Eq. (lA24p to proof the orthonormality (|3T|) of the Weyl operators (l30l) 
TrUl^Uij = ^ ^ e^^^'-'^") (p|(A; + m)modc/)(A;|fc)((fc + j)modc/|p) 

P=0 fc,A;=0 

= ^ J] (p| (A; + m) mod rf) ((A; + j) mod d\p) 6^-^ 

P=0 k,k=0 
k=0 



= d5ni5mj- (A25) 

4. Expansion into WOB 

Formula fl54l) for the Bell state in terms of WOB we derive in the following way. We 
express the standard matrices by the WOB fl371) . rewrite the indices and separate the non- 
vanishing terms 

j,k=i 

d-l d-1 

^ ^ Z-^ ^ Ul(k-j)modd ® Uv{k-j)modd 

j,k=Ol,l'=0 



Ci3 

m,fc=0 /,/'=0 



V m m l,l'-l+l'=d 

+ i E E f E ^'"^ f/^^ f/^,^ . (A26) 



The last term in Eq. (1A26|) vanishes due to relation (lA24p . Identifying Uqo = 1 and using 



the notation with negative values of the index /, which have to be considered as mod d , we 
gain the formula 

1 1 '^'^ 

k+>«l = ^1®! + ^ E^'™®^"'™' (/,"^) 7^ (0,0). (A27) 

Lm=0 
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